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Abstract — A deterministic attitude estimator for a rigid 
body under an attitude dependent potential is studied. This 
estimator requires only a single direction measurement to 
a known reference point at each measurement instant. The 
measurement cannot completely determine the attitude, but 
an attitude estimation scheme based on this measurement is 
developed; a feasible set compatible with the measurement is 
described and it is combined with an attitude dynamics model to 
obtain an attitude estimate. The attitude is globally represented 
by a rotation matrix, and the uncertainties are described by 
ellipsoidal sets. A numerical example for a spacecraft in a 
circular orbit is presented. 

I. Introduction 

The attitude of a rigid body is defined by the orientation 
of a body-fixed frame with respect to a reference frame, and 
the attitude is represented by a rotation matrix that is a 3 x 3 
orthogonal matrix with determinant of 1, which transforms 
a representation of a vector in a body-fixed frame into one 
represented in the reference frame. Rotation matrices have a 
group structure denoted by S0(3). In spacecraft applications, 
the attitude is usually determined by using a set of direction 
measurements. The directions to objects such as the sun, 
stars, and geomagnetic fields, assumed to be known in the 
reference frame, are measured in the body-fixed frame in 
order to determine the rotation matrix. 

Attitude determination using multiple direction measure- 
ments with least squares estimation is known as Wahba's 
problem [1]. The original solution of Wahba's problem 
is given in [2], and solutions are expressed in terms of 
quaternions [3], and in terms of a rotation matrix [4]. 
Based on these attitude determination schemes, attitude 
estimation problems are studied in [5], [6] and [7]. The 
attitude determination/estimation procedures using Wahba's 
problem formulation require at least two different direction 
measurements at each measurement instant. This places a 
stringent constraint on the dynamic estimation of spacecraft 
attitude. 

A single direction measurement provides some informa- 
tion about the attitude; it is guaranteed that the rotation 
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matrix lies in a one dimensional subgroup of the three 
dimensional special orthogonal group S0(3), which is dif- 
feomorphic to the one-sphere S^. The attitude is not com- 
pletely determined at a single measurement instant. If the 
process is coupled with an attitude dynamics model, an 
attitude estimation scheme can be developed using single 
direction measurements. An attitude determination scheme 
using single direction measurements is studied in [8], but 
this approach requires an additional arc length measurement. 

Most existing attitude estimation schemes use generalized 
coordinate representations of the attitude. As is well known, 
minimal coordinate representations of the rotation group, 
such as Euler angles, lead to singularities. Non-minimal 
coordinate representations, like the quaternions, have their 
own associated problems. Besides the extra constraint of 
unit norm that one needs to impose on the quaternion, the 
quaternion representation, which is diffeomorphic to SU(2), 
double covers S0(3). So, it has an inevitable ambiguity in 
expressing the attitude. 

A stochastic state estimator requires probabilistic models 
for the state uncertainty and the noise. However, statistical 
properties of the uncertainty and the noise are often not 
available. An alternative deterministic approach is to spec- 
ify bounds on the uncertainty and the measurement noise 
without an assumption on their distribution. Noise bounds 
are available in many cases, and deterministic estimation 
is robust to the noise distribution [9]. An efficient but 
flexible way to describe the bounds is using ellipsoidal 
sets, referred to as uncertainty ellipsoids. The deterministic 
estimation process is based on set theory results developed 
in [10]; optimal deterministic estimation problems using the 
uncertainty ellipsoids are studied in [11] and [12]. 

In this paper, a deterministic attitude estimator which 
requires a single direction measurement at each measure- 
ment instant is presented. A feasible set in S0(3) that 
is compatible with the measurement is represented by Lie 
algebra elements and the exponential map. It is compared 
with the attitude dynamics model to obtain an updated 
estimate. The estimation scheme presented in this paper has 
the following distinctive features: the estimator requires only 
a single direction measurement at each measurement instant, 
the attitude is represented by a rotation matrix without any 
local parameterization, and the deterministic estimator is 
distinguished from a Kalman or extended Kalman filter 

This paper is organized as follows. The attitude dynamics 
and uncertainty model are given in Section II. The attitude 
determination scheme and the attitude estimation scheme 
using single direction measurements are presented in Section 
III and IV, which is followed by a numerical example in 



Section V. 

II. Attitude dynamics and Uncertainty model 

A. Equations of motion 

We consider estimation of the attitude dynamics of a rigid 
body in the presence of an attitude dependent potential, 
[/(•) : S0(3) M, i? e S0(3). Systems that can 
be so modeled include a free rigid body, spacecraft on a 
circular orbit with gravity gradient effects [13], or a 3D 
pendulum [14]. The continuous equations of motion are 



Jfl + flxJfl = M, 
R = RS{n), 



(1) 
(2) 



where J e ffi^^^ is the moment of inertia matrix of the rigid 

body, f7 G is the angular velocity of the body expressed 
in the body-fixed frame, and S{-) : i-^ so(3) is a skew 
mapping defined hy S{x)y = a; x y for all a;,y e M^. 
The vector M G R"^ is the moment due to the potential, 
determined by S{M) — ^ R — R^^^ or more explicitly. 



M = ri X Vr^ + r2 X Vr^ + rs X Vr^^ , 



(3) 



where ri,Vri G M^^" are the ith row vectors of R and 
respectively. 

General numerical integration methods like the popular 
Runge-Kutta schemes, typically preserve neither first in- 
tegrals nor the characteristics of the configuration space, 
SO (3). In particular, the orthogonal structure of the rotation 
matrices is not preserved numerically. It is often proposed to 
parameterize (2) by Euler angles or quaternions instead of 
integrating (2) directly. However, Euler angles yield only lo- 
cal representations of the attitude and they have singularities. 
Unit quaternions do not exhibit singularities, but they have 
the manifold structure of the three-sphere and double 
cover SO (3). Consequently, the unit quaternion representing 
the attitude is inevitably ambiguous. In addition, general 
numerical integration methods do not preserve the unit length 
constraint. Therefore, quaternions have the same numerical 
drift problem as rotation matrices. 

Lie group variational integrators preserve the group struc- 
ture without the use of local charts, reprojection, or con- 
straints, they are symplectic and momentum preserving, and 
they exhibit good energy behavior for an exponentially long 
time period. The following Lie group variational integrator 
for the attitude dynamics of a rigid body is presented in [14]: 

hSiJrik + ^Mfe) =FkJd- JdF^, (4) 
^fe+i = RkFk, (5) 
Jfifc+i = Fj^mk + ^F^Mk + ^Mfc+i , (6) 

where G R"^ ^ ^ is a nonstandard moment of inertia matrix 
defined by Jd = itr[J] hy,z - J, and Fk G S0(3) is 
the relative attitude between integration steps. The constant 
h G M is the integration step size, and the subscript k 
denotes the /cth integration step. This integrator yields a map 
{Rk,flk) (-Rfe+i,^^fe+i) by solving (4) to obtain Fk G 



SO (3) and substituting it into (5) and (6) to obtain Rk+i and 
fik+i- The only impUcit part is (4). The actual computation 
of Fk is done in the Lie algebra so(3) of dimension 3, and 
the rotation matrices are updated by multiplication. So this 
approach is distinguished from integration of the kinematics 
equation (2), and there is no excessive computational burden. 
We use these discrete equations of motion to propagate 
the attitude dynamics between measurements during the 
estimation process. 

B. Uncertainty Ellipsoid 

We describe uncertainties of the attitude dynamics by 
using ellipsoidal sets referred to as uncertainty elUpsoids. 
An uncertainty ellipsoid in R" is defined as 



^1 



{x,P) = |.T G M" (x - x fp-\x - x) < l| 



(7) 



£{R,n,P) = \ R€S0{3),nG 



where x G M", and P G R"^" is a symmetric positive 
definite matrix. We call x the center of the uncertainty 
ellipsoid, and P is the uncertainty matrix that determines 
the size and the shape of the uncertainty ellipsoid. The size 
of an uncertainty ellipsoid is measured by tr[P] which is the 
sum of the squares of the semi principal axes of the ellipsoid. 

The attitude dynamics evolves on the 6 dimensional tan- 
gent bundle, TS0(3). We identify TS0(3) with SO (3) x 
so (3) by left trivialization, and we identify so (3) with R"^ 
by the isomorphism S{-). The uncertainty ellipsoid centered 
at {R, O) G TS0(3) is induced from an uncertainty ellipsoid 
in R*^; 

e£^e{Oe,P)Y 
(8) 

where S{C) = logm (^FfR^ G so(3), Sn = Q-Q gM.^, 
and P G R®^® is a symmetric positive definite matrix. An 
element {R, f2) G £{R, fl, P) can be written as 

P = Pexp5(C), n = Q + Sn, 

for some x = [C; SQ] G R® satisfying x^ P~^x < 1. 

We assume that the initial conditions are bounded by a 
prescribed uncertainty ellipsoid 

(Po,r!o) Gf(i?o,f^o,Po), (9) 

where Pq g R^^® is a symmetric positive definite matrix that 
defines the shape and the size of the uncertainty ellipsoid. 

III. Attitude Determination with a Single 

Direction Measurement 

In the attitude determination problem, we measure direc- 
tions to points in the reference frame. We assume that the 
directions to these points are known in the reference frame. 
This either requires that the points are located far away from 
the spacecraft or the relative location of the spacecraft is 
known exactly. The directional sensor is fixed in the body- 
fixed frame, and the measurements are representations of the 
direction vectors in the body-fixed frame. The representations 
in the body-fixed frame are transformed into those in the 
reference frame by multiplication with the rotation matrix 
that defines the attitude of the rigid body. 



A. Exact measurement 

Let the direction to a known point in the reference frame 
be e G S^, and let the corresponding vector represented in 
the body-fixed frame be 6 G We first assume that the 
direction measurement has no error, so the direction b is 
exact. Since we only measure a direction to a point, we 
normalize e and b so that they have unit lengths. The vectors 
e and b are different representations of the same vector from 
the spacecraft to the known point, and they are related by 
a rotation matrix R e SO (3) that defines the attitude of the 
rigid body 



e = Rb. 



(10) 



This equation provides a two-dimensional constraint on 
the three-dimensional rotation matrix. Consequently, a single 
direction measurement does not completely determine the 
attitude. This corresponds to the fact that if we rotate the 
rigid body about the direction e in the reference frame, then 
the measured direction b is not changed. The rotation matrix 
has one-dimensional uncertainty represented by any rotation 
about the direction e in the reference frame, or equivalently, 
any rotation about the direction b in the body-fixed frame. 

Suppose that R° G SO (3) is a particular rotation matrix 
satisfying (10). This rotation matrix can be represented in 
several ways. For example, if b and e are not co-Unear, 

bx e 



R° = exp 



cos-^(6^e) S 



exp 5(6*° 6), (11) 



||6 X e|| 

where the constant 9° E can be arbitrarily chosen. The 
rotation matrix that represents the attitude of the rigid body 
can be written in terms of R° as 



R = R° exp [9S{b)] 



(12) 



for a 61 G 

In summary, if the single direction to a known point is 
measured exactly, the rotation matrix lies in the following 
one dimensional subgroup of S0(3): 



Rg{R° exp [eS{b)] 



6* G §^ 



(13) 



B. Measurement error 

We now consider the effects of small measurement errors. 
Let 6 G §^ be the measured direction of the direction b. 
Since we only measure directions, we normalize b and h so 
that they have unit lengths. Therefore it is inappropriate to 
express the measurement error by a vector difference. The 
measurement error is modeled by rotation of the measured 
direction 



b = exp [S{v)] b, 
^b + S{v)b, 



(14) 



where G M'^ is the Euler axis of rotation from b to b, and 
is the corresponding rotation angle error in radians. The 
measurement error is bounded by an uncertainty ellipsoid 



z/g£m3(03,5) 



for a symmetric positive definite matrix S G M?'^^. The 
magnitude of the measurement error is assumed to be small. 

Let -R° G S0(3) be a rotation matrix obtained by (11) for 
the measured direction b. We express the difference between 
R° and R° using the exponential map: 



R° = ^°exp[5(C°)] 



(16) 



for some C° G Since we make the small measurement 
error assumption, the norm of the vector C,° is considered to 
be much smaller than tt, i.e. ||C°|| <C tt. Since e = R°b = 
R°b, we obtain 

R°b = R°b, 

= R° exp[5(C°)] {hy3 + S{v)}b, 

Thus we have S{(° + v)b = 0, which is equivalent to 

(l°=cb-v (17) 

for any constant c G M. Since ||C°|| < |c| -|- < tt, the 
constant c is smaller than tt, i.e. |c| -C tt. 

Substituting (14), (16), and (17) into (12), we obtain 



R = R° exp S{cb - v) exp 9S ((/sxs + S{v))b) 



(18) 



(15) 



for constants c and ^ G S^. 

In summary, if the single direction measurement has a 
small error represented by (14), then the attitude of the rigid 
body can be written in terms of the measured direction b and 
the measurement error v as (18). This expression includes 
the uncertainty caused by the measurement error as well 
as the uncertainty due to the single direction measurement 
assumption. The constant 9° G to determine R° and R° 
is specified by the following estimation procedure. 

IV. Attitude Estimation with a Single Direction 
Measurement 

The deterministic estimation scheme using uncertainty 
ellipsoids is introduced first. A deterministic estimator for the 
attitude and the angular velocity of a rigid body is developed 
by using the preceding attitude determination scheme. 

The subscript k denotes the fcth discrete index. The su- 
perscript / denotes the variables related to the flow update, 
and the superscript m denotes the variables related to the 
measurement update. ~ denotes a measured variable, and • 
denotes an estimated variable. 

A. Deterministic estimation 

We use deterministic bounded estimation using ellipsoidal 
sets, referred to as uncertainty ellipsoids, to describe the un- 
certainty and measurement noise. The estimation process has 
three steps similar to those in the Kahnan filter: prediction, 
measurement, and filtering steps. We assume that the initial 
condition lies in a prescribed uncertainty ellipsoid, which 
is propagated in time using the equations of motion. This 
defines a prediction step. The measurement error bound is 
described by a measurement uncertainty ellipsoid. Then we 
can guarantee that the state lies in the intersection of the 




(a) Propagation of uncertainty ellipsoid (b) Filtering proce- 

dure 

Fig. 1. Uncertainty ellipsoids 

predicted uncertainty ellipsoid and the measured uncertainty 
ellipsoid. The intersection of the two ellipsoids is an irregular 
shape, which is not efficient to compute and store. Instead 
we find a minimal ellipsoid that contains this intersection. 
This procedure is repeated whenever new measurements are 
available. 

This deterministic estimation procedure is illustrated in 
Fig. 1 . The left figure shows time evolution of an uncertainty 
ellipsoid, and the right figure shows a cross section at a 
fixed measurement instant. At the fcth time step, the state 
is bounded by an uncertainty ellipsoid centered at Xk- This 
initial ellipsoid is propagated through time. Suppose that the 
state is measured next at the {k + l)th time step. In this 
single direction measurement estimation, the measurement 
ellipsoid degenerates to a strip. At this instant, the actual 
state lies in the intersection. In the estimation process, we 
find a new ellipsoid that contains this intersection, as shown 
in the right figure. The center of the new ellipsoid, Xk+i 
gives a point estimate of the state at time step {k + I), and 
the magnitude of the new uncertainty ellipsoid measures the 
estimation accuracy. The deterministic estimates are optimal 
in the sense that the sizes of the ellipsoids are minimized. 

B. Flow update 

Suppose that the attitude and the angular momentum at 
the fcth step lie in a given uncertainty ellipsoid; 

{Rk,nk) e£{Rk,nk,Pk), 

and a new measurement is taken at the (fc + ^)th time step. 

The flow update finds the center and the uncertainty 
matrix that define the uncertainty ellipsoid at the (fc + ^)th 
step using the given uncertainty ellipsoid at the fcth step. 
Since the attitude dynamics of a rigid body is nonlinear, the 
admissible boundary of the state at the (fc + ^)th step is not 
an ellipsoid in general. We assume that the given uncertainty 
ellipsoid at the fcth step is sufficiently small that attitudes 
and angular velocities in the uncertainty ellipsoids can be 
approximated using the linearized equations of motion. Then 
we can guarantee that the uncertainty set at the (fc + ^)th 
step is an ellipsoid, and we can compute its center and its 
uncertainty matrix at the (fc + /)th step separately. 



Center: For the given center at step fc, [Rk, ^k), the center 
of the uncertainty ellipsoid at step (fc + /) is (^^^^£7^^;) 
obtained using the discrete equations of motion, (4), (5), and 
(6): 

hSiJdk + ^Mk) = FkJd - JdF^. (19) 
= Rkh. (20) 
mi^^ = Fk^^lk + ^Fk^Mk + ^Mk+i- (21) 

This integrator yields a map {Rk, ^k) '-^ i^l+i ' ^fc+i)' 
this process is repeatedly applied to find the center at the 
(fc + Oth step, {Ri^i,Cli^i). 

Uncertainty matrix: We assume that an uncertainty el- 
lipsoid contains small perturbations from the center of the 
uncertainty ellipsoid. Then the uncertainty matrix is propa- 
gated by using the linearized flow of the discrete equations 
of motion. At the (fc + l)th step, the uncertainty ellipsoid is 
represented by perturbations from the center iRl_^_i,^^l_^_i) 
as 

Rk+i^RU,expS{Cl+,), 

for some Cl+i'^^l+i ^ The uncertainty matrix at 
the (fc + l)th step is obtained by finding a bound on 
d+iT^^i+i £ Assume that the uncertainty ellipsoid 
at the fcth step is sufficiently small. Then, Cfc+ii'^^fc+i 
are represented by the following linear equations using the 
results presented in [13] 

where Xk = [Ck^Snk] £ R^, and A| e R^xe ^.^^ suitably 
defined. Since {Rk,^k) e £{Rk,^k,Pk), Xk G 5]r6(0, Pfc) 
by the definition of the uncertainty ellipsoid given in (8). This 
implies that A^Xk lies in the following uncertainty ellipsoid 




Thus, the uncertainty matrix at the (fc + l)th step is given 
by 

<i-4^^(4)'^- (22) 

The above equation is then applied repeatedly to find the 
uncertainty matrix at the (fc + /)th step. 

In summary, the uncertainty ellipsoid at the (fc + Z)th step 
is computed using (19), (20), (21), and (22) as: 

{Rk+i,nk+i)e£{Ri^^,ni^i,Pi^^), (23) 

C. Measurement update 

The measurement update finds an uncertainty ellipsoid in 
the state space using the measurement and the measurement 
error models described in Section III. A feasible set of 
rotation matrices that is compatible with the single direction 
measurement is described in (18). We find an expression for 



the measurement uncertainty ellipsoid such that it contains 
the set described by (18). 

Elements in the measurement uncertainty ellipsoid are 
expressed as 



(24) 



for the center i?™^, e S0(3) and some C^, e M^. We 
omit the subscript (fc + 1) hereafter for convenience, and it is 
assumed that the direction is measured at the (A: + l)th step. 

Center: Comparing (18) and (24), we choose the center 
of the measurement uncertainty ellipsoid as 



R"' = R°, 



exp 



cos-^(Fe) si 



b X e 
116 X ell 



expS{e°b), (25) 



for the constant 6° € which is determined by the 
following filtering procedure. 

Uncertainty Matrix: From (18) and (24), we have 

e^pSiC) exp[-5(c6 - u)\ = exp[^5 ({hxs + S{u))b) 

Since the vectors and C° = cb — u are assumed to be 
small, the above equation is approximated as 

C" - {cb -v)=eb + eS{u)b, 

which can be rewritten as 

- (61 + c)b = -v- eS{b)v, 

= {e- l)v - ^(73x3 + S{b))y. 

Since v S £r3(03, S) and ^ G S^, the terms in the right hand 
side satisfy 

where A"' = /sxa + S(b) e R^^^^. Therefore, the vector 
— + c)6 lies in an ellipsoid containing the vector sum 
of the above two ellipsoids. The expressions for the minimal 
ellipsoid containing the vector sum of two ellipsoids are 
given in [11]. Using the results, we have 



C-(0 + c)6efM3(O3,Po'"), 



(26) 



where 

Pi 



-{l + q-^)Q^ + {l + q)Q^, q=^ 

Qi = (1 + Tr)^^, = tt'^A"'-^SA 



From (26), we can guarantee that the vector i^"* lies in an 
ellipsoid containing the following union of the sets 

£R3(-7rfe,Po")U£R3(7r6,Po"). 

This is a consequence of the fact that an elUpsoid is convex 
and the assumption |c| -C tt. The ellipsoid that contains the 



union of two ellipsoids is obtained numerically by the LMI 
approach presented in [15]. 

fM3 (O3, P"*) D (fM3 (-7r6, Po"*) U {-Kb, Po"*)) . (27) 

In summary, a single direction measurement with small 
error guarantees that the rotation matrix is expressed as (24), 
where the center is given by (25), and the vector lies 
in the uncertainty ellipsoid given by (27). 

D. Filtering procedure 

The filtering procedure finds a new uncertainty ellipsoid 
compatible with both the predicted uncertainty ellipsoid and 
the measured uncertainty ellipsoid. The intersection of two 
ellipsoids is generally not an ellipsoid. We find a minimal 
uncertainty ellipsoid containing the intersection. 

The predicted uncertainty ellipsoid is based on W and 
the measurement ellipsoid is based on P™. In the following 
development, we assume that the difference between the ro- 
tation matrices W and P™ is small. Here we find a value of 
9° e at (25) such that the difference is minimized. Define 
an index J = Xi J3X3 — Rf^'^R"^ . A standard variational 
approach with the use of Rodriguez formula shows that the 
index is minimized when 



tr 



■ tan 



tr P/.TpAS'(6)2 



tx\R^''^R''S{b)\ sinr - tr 



Rf'^R^'Sib)' 



cose° > 0, 



where R^ e S0(3) is the first exponential of (25). The first 
equation is obtained by the optimality condition = 0, 
and the second inequality is obtained by > 0. These 

conditions define the value of 6^ G uniquely. 

We find a minimal ellipsoid containing the intersection 
of the predicted uncertainty ellipsoid and the measurement 
uncertainty ellipsoid. An element in the predicted uncertainty 
ellipsoid, {Rf,Clf) G £{R^ ,df ,pf), can be written as 

Rf = RfeS(C'), (28) 

Qf = nf + sn^, (29) 

for some iC^'jnf) e £se.{OexuP^)- We find an equivalent 
expression based on the measurement elUpsoid center P™. 
Define C"-'' G such that 

Rf = R"'e^''^""l (30) 

Thus, ("^^ represents the difference between the centers of 
the two ellipsoids. Substituting (30) into (28), 

Rf =RmeS(C-^)eS{C')^ 

~P™e^(^""''+f'), 

where we assumed that CJ^^ are sufficiently small to 
obtain the second equality. Thus, the uncertainty ellipsoid 
obtained by the flow update, £{W ,P-f) is identified by 
the center (P™, Cl^) and the uncertainty ellipsoid in R^. 

{rf,snf)e£^.{x^f,p^), 



where f™/ = [C^^jOa], and 5(C'"^) = logm (^BJ^'^Rf^ e 

so(3), snf = nf -nf e r^. 

We seek a minimal ellipsoid that contains the intersection 
of the following uncertainty elhpsoids. 

(£:R3(03,P™)f|^K^(^'"^^^)) C£Re{x,P), (31) 



Expressions for a minimal 



where x 

ellipsoid containing the intersection of two ellipsoids are 
presented in [11]. Using those results, x and P are given 
by 

X = (/6X6 - LH)x"'f, 

P = p{r)[{I - LH)pf{I - LHf + r-^LP'^L^], 
where L e R^^^, H e R?^^, and /3(r) e R are given by 
L = pfH'^[HPfH^ + r-ip'"]-\ 

H = [-^3x3, 03x3]^, 

/3(r) = 1 + r - {x"'f)'^H'^[HpfH'^ + r-^P"']-^Hx"'f , 

for a constant r, which is chosen such that tr[P] is minimized. 

In summary, a new uncertainty ellipsoid at the {k + /)th 
step is given by 



(-Rfe+lj^fe+l) € f (iife+1, Pfe+i), 



where 



P, 



fc+i 



P 



= + 5(1, 
Pfe+i = P. 



(32) 

(33) 
(34) 
(35) 



The entire procedure is repeated whenever a new measure- 
ment is available. 

The steps outhned above define a dynamic filter. The 
center of the uncertainty ellipsoid is considered as a point es- 
timate of the attitude and the angular velocity at the (fc + Z)th 
step. The uncertainty matrix represents the characteristics 
of the uncertainty, and the size of the uncertainty matrix 
represents the accuracy of the estimate. If the size of the 
uncertainty ellipsoid is small, we conclude that the estimate 
is accurate. This estimation is optimal in the sense that the 
size of the filtered uncertainty elhpsoid is minimized. 

E. Properties of the estimator 

The notable feature of this attitude estimator is that it 
requires a single direction measurement. Current attitude 
estimators based on the solution of Wahba's problem require 
at least two direction measurements at each instant. A single 
direction measurement provides only a two-dimensional con- 
straint for the six-dimensional tangent bundle. The informa- 
tion obtained from the attitude dynamics is utilized, together 
with the measurement, in order to estimate the attitude and 
the angular velocity of the rigid body. In this paper, it is 
assumed that the angular velocity is not measured, but the 
current results can be readily extended to incorporate angular 
velocity measurements. 



This attitude estimator has no singularities since the atti- 
tude is represented by a rotation matrix, and the geometric 
structure of the rotation matrix is preserved since it is updated 
by the structure-preserving Lie group variational integrator. 
The presented estimator can be used for highly nonhnear 
large angle maneuvers of a rigid body. It is also robust 
to the distribution of the measurement noise since we only 
use ellipsoidal bounds on the noise. The measurements need 
not be periodic, the estimation is repeated whenever new 
measurements become available. 

V. Numerical Example 

Numerical simulation results are given for the estimation 
of the attitude dynamics of an uncontrolled rigid spacecraft in 
a circular orbit about a large central body, including gravity 
gradient effects. The details of the on orbit spacecraft model 
are presented in [13]. 

The mass, length and time dimensions are normaUzed 
by the mass of the spacecraft, the maximum length of the 
spacecraft, and the orbital angular velocity, respectively. The 
moment of inertia of the spacecraft is chosen as J = 
diag[l, 2.8, 2]. The maneuver is a large attitude change 
completed in a quarter of the orbit. The initial conditions 
are chosen as 



Po — 

Po = -^3x3 



0.707 -0.707 
0.707 0.707 
1 



^0 = [2.32, 0.45, -0.59], 

do = [2.12, 0.55, -0.89]. 



The corresponding initial estimation errors are ||Co|| = 
45 deg, ||(5rio|| = 21.43 rad/s. Note that the actual initial 
attitude is opposite to the estimated initial attitude. The initial 
uncertainty matrix is given by 

Po = diag [2.28, 2.28, 2.28, 0.82, 0.82, 0.82] , 

so that x^Pq^xo = 0.857 < 1. 

We assume that the measurements are available twenty 
times. The inertial direction e to a known point is chosen 
from the columns of the following matrix. 



E = 



1 
1 




0.7071 

0.7071 

1 



-0.7071 
0.7071 




0.5 
0.5 
0.7071 



-0.5 
0.5 
0.7071 



A simple adaptive scheme is developed to choose the best 
inertial direction as the spacecraft rotates. The uncertainty 
matrix for the measurement noise is given by 



TT \ ' 

0.2 

180/ 



The direction measurement noise is normally distributed in 
the simulation. 

Fig. 2 shows simulation results, where the left figure 
shows the attitude estimation error and the angular velocity 
estimation error, and the right figure shows the size of the 
uncertainty ellipsoid. The estimation errors and the size 
of the uncertainty are reduced rapidly after the first few 
measurements; the estimation error for the angular velocity 
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= 0.2 
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(a) Estimation error 
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(b) Size of uncertainty tr[Pj. 
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Fig. 2. Estimation en'ors and uncertainty 



converges relatively slowly since the angular velocity is 
not measured directly. The terminal attitude error, and the 
terminal angular velocity error are less than 2.3 deg, and 
0.08rad/s, respectively. 

VI. Conclusions 

A deterministic attitude estimator for a rigid body under 
an attitude dependent potential is developed. This estimator 
requires only a single direction measurement to a known 
reference point at each measurement instant. A feasible set 
of rotation matrices compatible with the measurement is 
described in terms of Lie algebra elements, and it is com- 
pared with an uncertainty ellipsoid obtained from an attitude 
dynamics model, in order to obtain an updated attitude 
estimate. The attitude is globally represented by a rotation 
matrix, and the geometric structure of the rotation matrix is 
preserved by using a Lie group variational integrator. 
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